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Abstract 

The calculation of perturbative corrections to the spectral moments ob- 
servable in hadronic r decays is reconsidered. The exact order-a^ results 
and the resummation procedure of Le Diberder and Pich are compared with 
a partial resummation of the perturbative series based on the analysis of 
so-called renormalon chains. The perturbative analysis is complemented 
by a mo del- independent description of power corrections. For the con- 
tributions of dimension four and six in the OPE, it is demonstrated how 
infrared renormalon ambiguities in the definition of perturbation theory 
can be absorbed by a redefinition of nonperturbative parameters. We find 
that previous determinations of QCD parameters from a measurement of 
spectral moments in r decays have underestimated the theoretical uncer- 
tainties. Given the present understanding of the asymptotic behaviour of 
perturbation theory, the running coupling constant can be measured at 
best with a theoretical uncertainty 5as(m'^) ~ 0.05, and the gluon conden- 
sate with an uncertainty of order its magnitude. Two weighted integrals of 
the hadronic spectral function are constructed, which can be used to test 
the absence of dimension-two operators and to measure directly the gluon 
condensate. 
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1 Introduction 



The r lepton is the only known lepton heavy enough to decay into hadrons. Its 
decays provide a unique environment to study hadronic weak interactions at low 
energies. Because of its inclusive character, the total r hadronic width is expected 
to be calculable in QCD using analyticity and the operator product expansion 
(OPE) Since nonperturbative contributions to the total width turn out 

to be strongly suppressed, a measurement of the r width is considered an excellent 
way to extract the strong coupling constant Us with high precision at low energies. 
It has also been pointed out that hadronic r decays can more generally be used 
to test QCD and to extract some of its nonperturbative parameters |p-|T^. The 
idea is to measure spectral moments built from weighted integrals of the invariant 
hadronic mass distribution. As the total width, such moments are calculable in 
QCD, but they are more sensitive to nonperturbative contributions. An analysis 
of the spectral moments can also be used to search for non-standard effects, such 
as power corrections of order (A/mT-)^ [|13|, which are absent in the standard 



OPE approach of Shifman, Vainshtein and Zakharov (SVZ) ||15|| . 

A QCD-based analysis of power-suppressed effects by measuring spectral mo- 



ments in hadronic r decays has been pursued by two experimental groups [|T^, p!7 
There is a theoretical obstacle to any such analysis, however. The separation of 
perturbative and nonperturbative effects in QCD is intrinsically ambiguous. In 
order to define and extract power corrections in a meaningful way, one has to 
control the perturbative contributions to sufficiently high accuracy. But it is 
known that QCD perturbation theory provides an asymptotic expansion, which 
is factorially divergent at large orders. Related to this behaviour are the so-called 
renormalon singularities in the Borel transform of a perturbative series with re- 
spect to 1/as |18[-|^. The definition of the (resummed) perturbative expansion 



by itself is ambiguous; only the sum of perturbative and nonperturbative contri- 
butions in the OPE is well defined [^. Recently, the study of renormalons has 
received renewed attention. Efficient techniques have been developed to resum 
so-called renormalon chains, i.e. terms of order /5o "^a" (where /3o denotes the first 
coefficient of the /9-function), to all orders in perturbation theory [^-[Q. These 
investigations are interesting, since they provide an estimate of the importance 
of higher-order terms in a series, and moreover they elucidate the structure of 
nonperturbative contributions, which have to be included in the OPE in order 
to cancel the ambiguities of resummed perturbation theory. They are also useful 
for estimating the uncertainty in finite-order perturbative calculations. 

In the present paper, we reconsider the theoretical analysis of hadronic r de- 
cays in the light of these developments. This extends some recent analyses of 
renormalon contributions to the r hadronic width In Sect. ^ we define 

the spectral moments and discuss existing perturbative predictions for them. In 

^The applicability of the OPE in r decays has, however, been questioned recently 0. 
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Sects. 1^ and ^ we discuss the resummation of renormalon contributions for eu- 
clidean current correlators and for the spectral moments. Section |^ is devoted to 
a model-independent analysis of power corrections. In Sect. ^ we show in detail 
how the ambiguities of perturbation theory, which are due to the first two infrared 
(IR) renormalons, can be absorbed into a redefinition of some nonperturbative 
parameters. A discussion of our results and a study of the feasibility of extract- 
ing the running coupling constant and some nonperturbative parameters from 
data are presented in Sect. |^. There we construct two weighted integrals of the 
hadronic spectral function, which can be used to test the absence of dimension- 
two operators in QCD and to measure the gluon condensate. Technical details 
of our calculations are relegated to three appendices. 



2 Spectral moments 

The theoretical description of inclusive hadronic r decays involves two-point cor- 
relation functions of flavour- changing vector and axial vector currents, V-^j = 
Qi'j^qj and A^- = qil^l^qj, where i and j are flavour labels {i = u, j = d or s). 
These correlators admit the Lorentz decomposition [T = V ot A) 

I J d'x e'^- (0| T {ri;.(x), rr,.(0)t} |0) = (g^'g^ - g^^q') n« + q'q"" ^^viq') ■ 

(1) 

The superscript (J) denotes the angular momentum in the rest frame of the 
hadronic final state. The total hadronic r decay rate, normalized to the electronic 
one, can be written as 

9 

mi;. 

^ _ T{t ^ Ur + hadrons) _ dRrjs) 

where 



ds m'i \ rn^ 



1 + — Im n(°+^) (s -it) Im n(°) (s - it 

mi J mi 

(3) 

is the inclusive hadronic spectrum, and 

n(^)(^) = ^ (I K.pna^+A(^) + I K/nl^V+A(^)) • (4) 

The factor 5'ew — 1.0194 accounts for electroweak radiative corrections PT|. 

Le Diberder and Pich have pointed out the usefulness of considering weighted 
integrals of the spectrum di?^/ds |T^. Such integrals are linear combinations of 
the spectral moments 



M 



(J) _ A7r{k + 1) 



m 



2k+2 



Jdss^lmU^-^\s-it); k>0, (5) 
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which contain the dynamical information that can be extracted from hadronic r 
decays.0 The weighted integrals 



(A; + l)(A; + 3)(fc + 4) f\ / s \^ dRr{s) 
lik — / ds 



36S 



EW 



ds 



(6) 



which are normalized such that their perturbative expansion is Rk = l+a;s/7r+. . ., 
are given by 



Rk 



(fc + 3)(fc + 4) ^(0+1) _ (fc + l)(fc + 4) ^^(0+1) ^ (fc + l)(A: + 3) ^^(o+i) 



M 



fc+3 



(fc + l)(/c + 3)(fc + 4) ^(0)^ ^ 2(A: + l)(fc + 4) ^(o)^ _ (fc + l)(fc + 3) ^(o) 



3(A; + 2) 



In particular, we note that 



2MT + ^MT -M 



(0) 



,(0) 



Using the analyticity properties of the correlators in (|lD, the spectral moments 
can be written as contour integrals in the complex s-plane. One finds 



" 2m J s 

|s|=mj 



1 - 



mf 



fc+i 



where 



4vr2s4^n(-^)(s) 
ds 



(9) 



(10) 



are the logarithmic derivatives of the correlators, which are ultraviolet (UV) fi- 
nite. These functions are analytic in the complex s-plane, with discontinuities on 
the positive real axis. The contour integrals have two attractive features. First, 



to perform the integration along a circle with radius \s\ 



m: 



requires knowl- 



edge of the correlation functions only for large (complex) momenta. Second, the 



integrand in (|^) vanishes for s 



mz 



where the contour touches the physical 



cut. Therefore, unlike the spectral functions themselves, their moments are ex- 
pected to admit a theoretical description in QCD, which can be organized as an 
expansion in powers of {K/irirY^ 
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One of the purposes of this paper is to reconsider the calculation of the pertur- 
bative contribution in this expansion, i.e. the term with n = 0. This contribution 
can be calculated by setting the current quark masses of the light quarks m, d 



^Here we shall not pursue the idea to disentangle vector from axial vector, or strange from 
non-strange modes, although in principle this may be possible. 



fc+3 • 

(7) 
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and s to zero. In this limit the flavour- changing currents are conserved, so that 
the correlators with J = in ([l|) vanish. Moreover, chiral invariance implies that 
the spectral functions in the vector and axial vector channels are the same. It 
follows thatQ 

D^^:,'\-s) = Di-s) , =0. (11) 

The perturbative expansion of the function D{—s), 

D{-s) = l + ±KJ^^y, (12) 

n=l V TT y 

is known to order a^. For n > 3, the coefficients Kn depend on the renormaliza- 



tion scheme. In the MS scheme, one finds (with nj = 3 quark fiavours) ||32||-i36 



Ki = 1, 
299 

i^2 = —-9C(3)~ 1.63982, 

MS ^ 58057 _ 779 75 ^ 

^ 288 4 2 ^ ' ^ ^ 

with C(3) ~ 1.20206 and C(5) ~ 1.03693. The next coefficient, K^, has been 
estimated using the principle of minimal sensitivity |^ and the effective charge 
approach with the result that |^ 



Kest ^ 27.5 . (14) 
The perturbative contributions to the moments in (^ are 

A^gp+^rl = M, , = , (15) 

with 

M, = -L / - (1 - x'^') D{-xml) = 1 + £ d„ ('^MlY . (16) 
2vrz J X ~[ \ 71 / 

\x\=l 

The coefficients dn can be obtained using the renormalization-group equation 
(RGE) 

^ djj,'^ ^ ^^'^^ )/5f4/^ )] = -as(/x ) 2^ Pn y j (17) 



■^We set |y„ciP + iKtsP = 1, which is an excellent approximation, since experimentally 
Vub? < 1.6 X 10-5. 
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Table 1: Perturb ative contributions to the moments at order a^. 





Mo 


Ml 


M2 


M3 


M4 


M5 


0.26 


1.116 


1.101 


1.096 


1.094 


1.093 


1.092 


0.28 


1.128 


1.110 


1.104 


1.102 


1.100 


1.099 


0.30 


1.141 


1.119 


1.113 


1.110 


1.108 


1.107 


0.32 


1.154 


1.128 


1.121 


1.118 


1.116 


1.114 


0.34 


1.168 


1.138 


1.129 


1.125 


1.123 


1.122 


0.36 


1.183 


1.148 


1.138 


1.133 


1.131 


1.129 


0.38 


1.198 


1.158 


1.146 


1.141 


1.138 


1.136 



for the running coupling constant. For n/ = 3, the first three coefficients of the 
^-function are po = 9, pi = 64 and /3p = 3863/6 The result is p 



di = 1, 

d2 = K2 + -ttA^ ^ 1.63982 ' ^'^^ 



4(fc + l) k + l' 



' ' 2(A; + 1) V^" 8 y 8 V(A; + 1)2 

11.3792 10.125 



-10.2839 + 



k+1 (k + l) 



2 



+ 4(fcTl)r°^ +y^' + ^~^j + 32(fc + l)3 

,^MS irrnrr 46.238 94.810 68.344 
Kr - 155.955 - + — r + 



A;+ 1 (A; + 1)2 ' (A; + 1)3 ■ ' ' 

As mentioned above, the exact result for the coefficient is unknown. Nev- 
ertheless, it follows that ratios and differences of moments are known exactly 
to order af, since this unknown coefficient cancels out. In Table |l], we show the 
order-a^ predictions for the first six moments as a function of the value of ^^(m^) 
in the MS scheme. 

Le Diberder and Pich have argued that one should improve the perturba- 
tive prediction for the spectral moments by performing a partial resummation 
of higher-order terms p, |T^. Their observation was that the evolution of the 



running coupling constant along the integration contour in (|l^) generates higher- 
order corrections that are expected to be larger than the "genuine" higher-order 
corrections to the correlator D{—s). These corrections can be resummed by writ- 
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Table 2: Perturbative contributions to the moments obtained from the re- 
summation PROCEDURE OF Le DiBERDER AND PiCH (LP) [|l2| . 













MLP 




0.26 


1.110 


1.090 


1.087 


1.086 


1.085 


1.084 


0.28 


1.120 


1.095 


1.093 


1.091 


1.090 


1.089 


0.30 


1.130 


1.100 


1.098 


1.096 


1.095 


1.094 


0.32 


1.140 


1.105 


1.103 


1.101 


1.100 


1.099 


0.34 


1.150 


1.109 


1.107 


1.105 


1.104 


1.103 


0.36 


1.160 


1.112 


1.112 


1.109 


1.108 


1.107 


0.38 


1.169 


1.114 


1.116 


1.113 


1.111 


1.111 



mg 



n=l 



= ^ / !!£(i_.^«)(^i(z^)". (19) 

2m J X V TT / 

\x\=l 

The contour integrals over the running coupling constant can be performed nu- 
merically (see Appendix |^). The effect of this resummation is significant. In 
Table ^, we show the results obtained using the two-loop RGE for the running 
coupling constant. The differences with respect to fixed-order perturbation the- 
ory are of order a few per cent. We shall comment on the usefulness of this 
resummation procedure at the end of Sect. ^. 

In this paper we investigate another type of resummation, which deals with 
the so-called renormalon- chain contributions |18H|30|. These are terms of order 



/3o~^a" in the perturbative series for D{—s) and for the moments M^, where 
/3o = 11 — |njis the first coefficient of the /3-function. A graphical representation 
of a renormalon chain is shown in Fig. |I]. The motivation to resum these terms to 
all orders is the following: The coefficients in the perturbative series (|12D and ([16D , 
which are polynomials in the number of quark flavours nj, can be reorganized as 
polynomials in /3o, i.e. 

dn = 5n + S^n-'^ Po~' + ■ ■ ■ + Po + ^1°^ • (20) 

Since in QCD the value of /5o is large, it is conceivable that a reasonable approx- 
imation to these coefficients is provided by the first term, i.e. the one with the 
highest power of /Sq. Some formal arguments supporting this assertion can be 
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found in Ref. |]30[. Below, we shall refer to this approximation as the "large- /3o 
limit". In this limit, one finds 

Kl = 1 , 

297 

«2/5o = ^-9C(3)~ 1.55649, 

. 48924 513 , 
"^^^o^^SS ^C(3)^ 15.7112, 

165537 49329 ^, , 10935 , 
al = —^^ ^ C(3) - C(5) ^ 24.8320 , (21) 



which can be compared with the exact values of Ki, K2 and given in (|Tt 
as well as with the estimated value of given in (0). With the exception of 
K^^, the large-/3o limit seems to work quite well. We also quote the corresponding 
results for the coefficients in the expansion of the moments. They are 

Si = l, 

/II 1 \ 2 25 

'^2/5o = (y - C(3) + j /5o ^ 1.55649 + 



4(A; + l)y'^^ ■ k + 1 

^3/5„^ ^ I f - ^C(3) - ^ + (H - ic(3)) ^ + ^ I PI 

which can be compared with the exact values of di, d2 and d^^ given in (|T^). 



Figure 1: Renormalon-chain contribution to the correlator 
D{—q^). The shaded bubble represents a self-energy insertion on 
the gluon propagator. We note that bubble summation is not a 

gauge-invariant procedure in a NON-ABELIAN THEORY; THE FIGURE IS 

meant as an illustration only. 

In Table ^, we show the results obtained in the large-/3o limit for the pertur- 
bative contributions to the moments at order a^. For the difference between the 
"large-/3o results" and the exact ones, we find 

- 1 (^)%(^-15«3,4c(5,.^)(H.)%0,„;). 

(23) 
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Table 3: Large-/5o approximation to the perturbative contributions to the 

MOMENTS AT ORDER a^. 





^large-/3o 




^^arge-/3o 


^large-/3o 




j^^arge-/3o 


0.26 


1.118 


1.104 


1.100 


1.098 


1.097 


1.096 


0.28 


1.131 


1.114 


1.109 


1.107 


1.106 


1.105 


0.30 


1.144 


1.124 


1.119 


1.116 


1.115 


1.114 


0.32 


1.158 


1.135 


1.128 


1.125 


1.124 


1.122 


0.34 


1.173 


1.146 


1.138 


1.135 


1.133 


1.131 


0.36 


1.189 


1.157 


1.149 


1.145 


1.142 


1.141 


0.38 


1.205 


1.169 


1.159 


1.155 


1.152 


1.150 



where as = Q;s(m^) in the MS scheme. A comparison of Tables ^ and ^ shows 
that these differences are rather small, supporting to some extent the usefulness 
of the large- /5o limit. 

In Refs. [^-[0, efficient techniques have been developed to resum the terms 
of order /?o~^'^s exactly to all orders in perturbation theory. Note that such a 
resummation includes the main part of the resummation of Le Diberder and Pich. 
In addition, a partial resummation of the higher-order coefficients Kn in the series 
( p!2D is achieved. This is important, since these coefficients are known to diverge 
factorially for large n, reflecting the asymptotic behaviour of perturbative expan- 
sions in QCD. The perturbative series for D{—s) and for the moments Mk are not 
Borel summable, implying that the definition of these quantities in perturbation 
theory has an intrinsic ambiguity, which has to be cured by adding nonpertur- 
bative contributions |24|. The resummation of renormalon-chain contributions 
allows one to quantify (or at least estimate) this ambiguity. Moreover, unlike any 
truncated perturbative series, the resummed expressions for D{—s) and for the 
moments are independent of the renormalization scheme. Nevertheless, one 
has to keep in mind that this method provides only a partial resummation of a 
perturbative series, and it is not obvious that it gives a better approximation to 
the full series than fixed-order perturbation theory. Therefore, to be conservative 
we shall take the differences between our procedure and others as an indication 
of the perturbative uncertainty. 



3 Renormalon resummation for the D function 

We start by considering the correlator D{—s) in the euclidean region, where 
—s = Q^> 0. The terms of order Pq'^o^ in the perturbative series ( |I2D can be 
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resummed by evaluating the Borel integral 



4 /■ , - _ / Atxu 



-DBorei(Q^) = 1 + ^ / duSniu) exp 



1 + - /dn5z,H(^) , (24) 







with the Borel transform PSI, EB] 



A is the scale parameter in the one-loop expression for the running coupling 
constant, and C is a scheme-dependent constant, with C = —5/3 in the MS 
scheme (C = in the V scheme, see below). The coefficients k„ have been 
defined in (pO|). The function Sd{u) contains pole singularities on the real u axis. 
The singularities on the negative axis are related to the UV behaviour of Feynman 
diagrams and are called UV renormalons. The singularities on the positive axis, 
which arise from the integration over low virtual momenta in Feynman diagrams. 



are called IR renormalons [|18| . In the present case there are UV renormalon poles 
located at m = —1, —2, . . . and IR renormalon poles at u = 2, 3, . . .. In general, a 
pole at u = Ui is associated with a factorial growth of the expansion coefficients 
Kn of the form r(n) (4Mj)~" (there is an extra factor n for a double pole), so 
that the renormalon singularity closest to the origin determines the asymptotic 
behaviour of the expansion coefficients. In the case of UV renormalons, these 
contributions have alternating sign and can be resummed by means of the Borel 



integral (|2^). However, for IR renormalons all terms have the same sign and the 



series is divergent. In fact, IR renormalon singularities fall on the integration 
contour in (0), making the value of the Borel integral ambiguous. A measure 
of the ambiguity is provided by the residue uf the nearest IR renormalon pole. 
In the present case this is located at m = 2 and leads to an ambiguity of order 
AVQ^ 

The appearance of IR renormalons indicates that perturbation theory is in- 
complete; it must be supplemented by nonperturbative corrections. Only the 
sum of all perturbative and nonperturbative contributions is unambiguous. For 
euclidean correlation functions of currents, the OPE provides a consistent frame- 
work for a systematic incorporation of nonperturbative effects In the case of 



the function Z)(Q^), such effects appear first at order l/Q^ and are parametrized 
by the gluon condensate |T3|, which has an ambiguity that compensates the am- 
biguity in the resummed perturbative series This will be discussed in 
more detail in Sect. |^. 

An efficient technique for evaluating the Borel integral (|23) has been developed 
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in Ref. where it was shown thatQ 



where 



-DBorel(Q^) 



1 A 

1 + - [-wnir)a{TQ'), 



Inr + ln(QVA 



(26) 



(27) 



Here ^^(rQ^) is the one-loop running couphng constant in the so-called V 
scheme, in which the coupling constant is defined in terms of the heavy-quark po- 
tential 1^^. The relation between the V scheme and the MS scheme is such that 
a^^\ix^) = as{e'-^ fi^) with C = —5/3. Thus, the parameter Ay can be obtained 
from the value of as(m^) in the MS scheme by 



Ay = rur exp (^^ 



2n 



(28) 



The function wd{t~) in (|26D describes the distribution of the virtuality of the 
gluon in the diagrams associated with the order-a^ correction to D{Q^). It is 
given by ||28[| 



32 

y 
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- Inr + r (1 



L2(-r) +lnr ln(l + r) 



r < 1 



32 [3 1 , , ^ 

— <^ - + - Inr + (1 + Inr) r 
3 l4 2 ^ ' 



+ r(l + r) 



L2(-r-i)-lnrln(l + r-^) 



r > 1 



(29) 



L2{x) = — /Q^yln(l — t) is the dilogarithm. In Fig. |, we show wd{t) as a 
function of Inr, which is the natural integration variable in (^). The expansion 
coefficients related to the moments of this distribution by 



dr 



r 



wn{T){-\nT-C) 



n-l 



(30) 



The value of the integral in (|2b| ) is independent of the renormalization scheme; 
changing the scheme simply amounts to rescaling the integration variable T.f\ 
However, as written above the integral is ambiguous because the integration 
contour runs over the Landau pole in the running coupling constant a(rQ^), which 
is located at r = rj;^ = Ay/Q^. This is nothing than the IR renormalon ambiguity 
mentioned above. Different regularization prescriptions for the Landau pole lead 



■*In Ref. HI], the function wd^t) was denoted twd^t). 
^This is true for regular (MS-like) schemes only. 
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to different results. A measure of the ambiguity is provided by the residue of the 
pole, i.e. 



The quantity ADj-^niQ'^) is equal to the sum of the residues of the IR renormalon 
poles in the integrand of the Borel integral in (^) Below we shall use 

the principle value prescription to regulate the Landau pole in (p6D, which is 
equivalent to calculating the principle value of the Borel integral. 

To illustrate the importance of the all-order resummation, we quote numer- 
ical results for the case = and asirn^) = 0.32 in the MS scheme. For 
these parameters, the order-a^ result in the large-/3o limit is D{m^) = 1.135 (the 
exact result is 1.126). The resummation of all terms of order /3o~^a" leads to a 
value that is significantly larger, DBorei("^r) — 1-151 ± 0.003. The central value 
corresponds to the principle value of the Borel integral, whereas the error is given 
by the renormalon ambiguity ADj-cn in (0)- It is instructive to investigate the 
behaviour of the perturbative series in the large-/?o limit. The expansion coef- 



ficients be obtained by expanding the Borel transform Sd{u) in (pSD in 

powers of u. For the MS and the V scheme, the results are given in Table ^. In 
the V scheme, the coefficients have alternating sign, and the nearest UV renor- 



malon at M = — 1 is dominant already at low orders |^0[. In the MS scheme, 
UV renormalons are not yet dominant at low orders due to an extra factor e^"''^ 
in the Borel transform, which suppresses the region of negative u. The first few 
coefficients have the same sign, and their growth is determined by the interplay 
of the UV renormalon at u = — 1 and the IR renormalon at u = 2. 

In Table ||, we show the partial sums D^^^ = 1 + J2n=i f^nPo'^ices/n)"' in 
the two schemes. Note that the asymptotic value, which is given by the integral 
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Table 4: Perturbative coefficients for the D function in the large-/5o limit. 



n 


1 


2 


3 


4 


5 


MS on-1 
«:„ Po 


1 


1.55649 


15.7112 


24.8320 


787.827 


'^n/^0 


1 


-2.19351 


18.1000 


-138.989 


1610.41 



in (|26[) , is scheme independent. In the MS scheme, the minimal term in the 
series is reached at n = 4 and gives a contribution of 0.3%. If one truncates 
the series at this point, the difference to the resummed result is 1.4%, which is 
about four times the renormalon ambiguity. In the V scheme, the minimal term is 
reached already at n = 2, and its contribution is 6.0%. If the series is truncated 
at this point, the difference to the resummed result is 4.6%. Since in the V 
scheme the leading UV renormalon is dominant, the error due to the truncation 
of the series is necessarily of order (Av/ttIt)^ — 6.7%. The fact that this error 
is much larger than the IR renormalon ambiguity, and moreover has a different 



power dependence on m,-, has been emphasized in Ref. |^3[. For comparison, we 
note that the leading nonperturbative contribution to the D function, which is 
proportional to the gluon condensate flSl , is given by (see Appendix P) 



I^power(m^) = ^ G^) ^ (1.2 ± 0.6)% . (32) 

Clearly, to give meaning to the value of nonperturbative parameters such as the 
vacuum condensates one has to control the higher-order behaviour of perturba- 
tion theory; in particular, it is necessary to resum the contribution of the first 
UV renormalon. If one believes that the large-/3o limit provides a reasonable 
description of the nature of this singularity, the resummation is achieved by per- 
forming the Borel integral, and then the residual ambiguity is due to the nearest 
IR renormalon and thus of order (Ay/mr)^. If, however, one wants to be more 
conservative and truncate the series at the last term that is exactly known, one 
is left with an uncertainty that is larger, namely of order the last term itself. 



4 Renormalon resummation for the moments 



Let us now turn to the resummation of renormalon- chain contributions for the 
spectral moments. In the Minkowski region, a simple representation of the Borel 
integral of the form (pBj) does not exist [|28[. Beneke et al. have shown that the 
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Table 5: Partial sums of the perturbative series for the D function in the 

LARGE-/3o LIMIT. In THE SECOND COLUMN WE SHOW FOR COMPARISON THE KNOWN 
EXACT RESULTS IN THE MS SCHEME. 







exact result (MS) 




1 


1.1019 


1.1019 


1.1648 


2 


1.1180 


1.1189 


1.1052 


3 


1.1346 


1.1256 


1.1863 


4 


1.1373 




1.0837 


5 


1.1459 




1.2795 


oo 


1.151 ±0.003 




1.151 ±0.003 



principle value of the Borel integral can be obtained from 27, 29 



Borcl 



1 °°dr 1 dr 

1±— / — W^fc(r) arctan[7ra(rm2)] + — Re [ —Wkir-ie), (33) 

nfJo J T (Jq J T 

-TL 



where tl = Ay/m^ is the position of the Landau pole in the coupling constant 
a{Tm1) defined in (|27|). The functions Wk{T) are related to weighted integrals of 
the distribution function Wd{t) |p8| : 



PVfe(r) = (A;± 1) jdxx^WD{r/x) 





(34) 



It is possible to perform these integrals explicitly. For /c = 0, we find that 

Inr 



39r 15 r 

W,{r) = ^ 4 - 3C(3) - ^ r - 2L,{-r) + [2r ± L^{-t) 



:i±r) 



L2(-r) ±lnr ln(l ± r) 



(35) 



32r 



1± — ±2L3(-r 
4r 



l-^-L,{-r-^)]\nr 



[l + r] 



L2{-r' ) - Inr ln(l ± r 



r > 1. 



L^{x) = Jq -J L2{t) is the trilogarithm. For k = 1 we obtain 

W,{t) = ^ { - 3C(3) - 2L3(-r)) ^ - - L,{-r) ) r Inr 
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+ (l + r) L2(-r) +lnr ln(l + r) 



r < 1 



(36) 



Mr \ 1 

W^{T) = —h + - + 2rL,{-r-' 



2 + ^ + rL2(-r-i)) Inr 
+ (1 + r) [i:2(-r-^) - Inr ln(l + t~^)\ ] 



r > 1. 



For the remaining cases, k > 2, we find 



A;3 



7fc2 - 18A; + 15 
^ 4(A;-1)3 ^ 

/I r (- 

+ - H ^ — 

\k k-1 k{k-l) 



4(A: + 1)2 
k-2 



4(A:-1)3 k{k-l) 



{k-iy 



r Inr 



k{k-l) 



L2(-r) +lnr ln(l + r) 
1 



T) + Sr>{-T) Inr - - InV 



r< 1, 



WJt) 



with 



32r 



(fc+1) 



P + A; + l , 3A; + 5 1 , /A; + 1 

1^ 



A;3 ^4(A; + l)2r 



+ 



2{k+ 1) r 



In- 



^\k^k-l ' k{k-l) 



L2i-T-')-\nT ln(l + 



-r 



1) 



r > 1 



(3^ 



These functions obey the normalization 



T 



(39) 



which yields the unit coefficient in front of tts/vr in the perturbative series for M^. 
Moreover, we note that Woo{t) = woij). In Fig. ^ we show Wkij) as a function 
of Inr for a few values of k. For k oo, these functions converge towards the 
function wd{t) depicted in Fig. ||. 

Weighted integrals of the functions H4(t) with powers of Inr determine the 
coefficients in the perturbative expansion of the moments Mk in the large-/5o limit. 
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1.4 




Figure 3: Functions Wkir) versus Inr, for A: = (solid line), k = 1 
(dashed line), k = 2 (dash-dotted line) and A; = 3 (dotted line). 
The asymptotic result for ^ oo is shown as the thin solid line. 



Expanding (|33D to order a^, one finds P8 



l+^-^((lnr)+C) 



+ ^((lnV)+2C (lnr)+C2-^ 



etc 



+0( 



TT 

(40) 



where = a<j(m^), C = —5/3 for tlie MS sclieme (C = for the V scheme), and 



(lnr) = -|-^lnrH^,(r)=4C(3) 



(InV) = - / — \n\WkiT) = 18 - 12C(3) + 

4 J T 



6 fc+ 1 
23 

y 



K(3) 



+ 



A; + 1 (A; + 1) 



#1) 



This does indeed reproduce the result given in (12^). 

In Table ^, we show the resummed results for the first six moments as a 
function of the value of the coupling constant as(m^) in the MS scheme, which 
according to (^) defines the scale parameter Ay. In order to perform the numer- 
ical integrations in (^3|), it is convenient to use the asymptotic behaviour of the 
functions Wkir) for large and small values of r, which is given in Appendix |B[ 
When expanded in powers of the coupling constant, our resummation reproduces 
the large-/?o limit for the perturbative coefficients given in (^2]). We can correct 
for the missing pieces, which are known up to order a^, using (|23|) . This leads to 
the numbers given in the lower portion of the table.Q 

^We note, however, that this correction depends on the renormalization scheme. The results 
presented here refer to the MS scheme. 
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Table 6: Predictions for the resummed perturbative contributions to the 

MOMENTS OBTAINED FROM THE PRINCIPLE VALUE OF THE BOREL INTEGRAL. ThE 
lower portion of the TABLE CONTAINS THE RESULTS CORRECTED FOR THE EXACT 
COEFFICIENTS UP TO ORDER af. 





^Bord 






^Bord 


^Borel 




0.26 


1.122 


1.097 


1.094 


1.092 


1.092 


1.091 


U.zo 


1 1 Qi^ 

i.iob 


T 1 no 
i.iUo 


i.iUU 


i.uyy 


i.Uyo 


i.uy 1 


0.30 


1.151 


1.108 


1.105 


1.104 


1.103 


1.103 


0.32 


1.168 


1.113 


1.108 


1.108 


1.108 


1.107 


0.34 


1.186 


1.117 


1.109 


1.110 


1.111 


1.111 


0.36 


1.205 


1.121 


1.109 


1.111 


1.112 


1.112 


0.38 


1.226 


1.126 


1.105 


1.109 


1.111 


1.112 


0.26 
0.28 
0.30 
0.32 
0.34 
0.36 
0.38 


1.119 
1.133 
1.148 
1.163 
1.180 
1.199 
1.218 


1.093 
1.098 
1.103 
1.106 
1.109 
1.112 
1.114 


1.090 
1.095 
1.098 
1.101 
1.100 
1.098 
1.093 


1.088 
1.093 
1.097 
1.100 
1.101 
1.100 
1.095 


1.087 
1.092 
1.097 
1.100 
1.101 
1.100 
1.097 


1.087 
1.092 
1.096 
1.099 
1.101 
1.101 
1.098 



The renormalon ambiguity in the perturbative contribution to the moments 
Mfc can be obtained by inspecting the corresponding Borel transforms 

Skiu) = -—- Soiu), 42 

k + 1 — U TTU 



or by performing the integrals |27 



1 ° d 

AMr = —lm [ —W,iT-te). (43) 



-TL 



The structure of IR renormalon poles is as follows: For k = Q there are single 
poles at M = 3, 4, . . .; for = 1 there are single poles at m = 2, 3, . . .; for = 2 
there is a double pole at m = 3 and singles poles at u = 4, 5, . . .; for > 2 there 
is a double pole ai u = k + 1 and single poles at u > 3 (m 7^ A; + 1). For the 
corresponding renormalon ambiguities, we obtain 



16 



Table 7: Perturbative coefficients for the first three moments in the large- 

/?0 LIMIT. 



n 


1 


2 


3 


4 


5 


k = 







1 


3.80649 


16.1854 


56.3140 


223.935 








1 


0.05649 


1.69928 


-17.9203 


139.694 


k = 


1 




1 


2.68149 


5.08959 


-35.7603 


-443.813 








1 


-1.06851 


-0.95907 


-32.6273 


154.154 


k = 


2 


rMS 


1 


2.30649 


2.51598 


-47.2767 


-424.721 






^iPr' 


1 


-1.44351 


-0.72018 


-31.0109 


207.944 



Let us again study the behaviour of the perturbative series in the large- /5o 
limit. The asymptotic growth of the expansion coefficients 5„ in (|20|), which can 
be obtained by expanding the Borel transform Sk{u) in powers of u, is determined 
by the UV renormalon at u = —1. However, the behaviour in low orders is 
governed by a complicated interplay of this UV renormalon with the nearest IR 
renormalon, which is located aX u = 2 for k = 1 and u = 3 for k ^ 1. For 
the first three moments, we show the coefficients 5„ in Table |^. In Table ^, we 
show the partial sums Mjf^^ = 1 + J2n=i ^nPo~^ («s/7r)" and compare them to 
the asymptotic value obtained from the resummation of renormalon chains. As 
before, the minimal term in the series is reached later in the MS scheme than in 
the V scheme. If one truncates the series at the minimal term, in both schemes 
the differences to the resummed results are typically of order 10~^ for k = 0, and 
of order 1-2% for k = 1,2. In all cases these differences are much larger than the 
renormalon ambiguities. 

An interesting question that can be addressed using the large-/5o limit is 
whether the resummation procedure of Le Diberder and Pich ||^ improves the 
convergence of the perturbative series. To answer this question, we investigate 
the partial sums 

N 

Mr'("^^ = l + E'^n/?o"-^/r^(m.), (45) 

n=l 

where the integrals I^^^{mT-) are evaluate using the one-loop /5-function for the 
running coupling constant (see Appendix This ensures that for ^ oo 
one recovers the large-/3o limit of the series. The partial sums for the first three 
moments are shown in Table ^. The results are interesting, as they indicate 
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Table 8: Partial sums of the perturbative series for the first three mo- 
ments IN the large-/3o limit. 





O.MS 




1,MS 




2, MS 




1 


1.1019 


1.1648 


1.1019 


1.1648 


1.1019 


1.1648 


2 


1.1414 


1.1663 


1.1297 


1.1358 


1.1258 


1.1256 


3 


1.1585 


1.1740 


1.1351 


1.1315 


1.1284 


1.1224 


4 


1.1645 


1.1607 


1.1312 


1.1074 


1.1234 


1.0995 


5 


1.1670 


1.1777 


1.1263 


1.1262 


1.1187 


1.1248 


oo 


1.1678 


± .00003 


1.1131 ±0.0039 


1.1080 ±0.0008 



that the convergence is not improved in a significant way. In some cases (such 
as A; = 0), fixed-order perturbation theory even converges better towards the 
asymptotic result. We believe that this observation is not specific for the large- 
Po limit, since the large vr^-terms resummed in the approach of Le Diberder and 
Pich are retained in this limit. 

Table 9: Partial sums of the perturbative series for the first three mo- 
ments obtained using the large-/?o limit of the resummation procedure of 
Le Diberder and Pich. 



N 


^LP^^) 
0,MS 




I.MS 




^LP^TV) 
2,MS 




1 


1.1155 


1.1852 


1.0973 


1.1350 


1.0940 


1.1316 


2 


1.1343 


1.1235 


1.1098 


1.1094 


1.1059 


1.1031 


3 


1.1526 


1.1903 


1.1180 


1.1104 


1.1141 


1.1208 


4 


1.1552 


1.1309 


1.1185 


1.1415 


1.1148 


1.1103 


5 


1.1622 


1.2018 


1.1175 


1.0503 


1.1155 


1.1341 


oo 


1.1678 


± .00003 


1.1131 ±0.0039 


1.1080 ±0.0008 



To summarize, in the first part of this paper we have investigated three 
schemes to calculate the perturbative contributions to the spectral moments M^,: 
fixed-order perturbation theory, the resummation procedure of Le Diberder and 
Pich, and resummed perturbation theory in the large-/3o limit. The results are 
summarized in Tables |I|, |^ and |[ As there is no strong argument to prefer one 
of these schemes over the others, the differences between the numerical results 
must be considered as theoretical uncertainties in the perturbative calculation of 
the moments. These differences are of order a few per cent. We note that taking 



ratios of moments (as proposed in Ref. |jT2[) does not improve the accuracy, al- 
though such ratios are known exactly to order af. For instance, the predictions 
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for the ratio Mi/Mq differ by a larger amount than the predictions for the indi- 
vidual moments. In the second part of the paper, we shall investigate what this 
uncertainty implies for the sensitivity to power corrections. 



5 Power corrections 

One of the goals of analysing moments of the r hadronic spectrum is to test QCD 
at the level of nonperturbative effects, which manifest themselves in the form 
of power corrections [jl2|. In the standard approach of SVZ, these corrections 
are parametrized in terms of vacuum expectation values of local, gauge-invariant 



operators, the so-called condensates |T^ . There have been attempts to extract the 
condensates of dimension four, six and even eight from the analysis of hadronic r 
decays 0-|n|, |]T^, [T^. However, the feasibility of such determinations is limited 



by the fact that, in order to be sensitive to power-suppressed effects, one has to 
control perturbation theory to a sufficient level of accuracy. 

There is also a conceptual need for an analysis of power corrections. As we 
have seen, the perturbative definition of the moments is ambiguous because of 
the presence of IR renormalons. In order to deal with these ambiguities in a 
consistent way, one is forced to add other nonperturbative contributions. Hence, 
we shall write 

M'-i!^ = Mi'l,, + A^g^r = M, 6j=, + . (46) 

Only the sum of the perturbative and nonperturbative contributions can be ex- 
pected to be well defined and unambiguous. We shall now present a model- 
independent analysis of power corrections, which does not rely on assumptions 
about their origin. Our goal is to relate the power corrections to the moments to 
those of the euclidean correlators D'^'^^Q'^) defined in (p!0D , for which we write 

D^'\Q') = dI'JAQ') + 4oL(Q^) = D{Q') 5j^, + t ^^1^ . (47) 

n=l ^ 

The nonperturbative quantities {0^2n{Q'^)) scale like A^". The only assumption 
we shall make is that these quantities have a weak, logarithmic dependence on 
. Writing then 

{0^2r^m) = exp (in^ TT^) (Off (m^)) = (%)^ {0^0) , (48) 

it is justified to treat V = d/dlnm^ as a small parameter. In fact, in QCD 
the nonperturbative quantities (Off (m^)) depend on rrir through the running 
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coupling constant and through the running quark masses,[] so that 



d d 
7^[a^(m^)] ml -^—^ - p[asiml)] as ^ 



^_^(2,„J_ + _„,_j. (49) 

Here 7m = — 2as/7r + ... is the anomalous dimension of the running quark mass. 
Inserting expression (^) into the contour integral (|^), we find 

OO -| 



-Mi^ower = E 4.n(V) {0^,i\ml)) , (50) 



n=l 



with 



1^1 = 1 



1 ^n+l sin(7rV) 

^ (/t + l-n + V)(n-V) vr ■ 



To second order in V, we obtain 



4,«(V) . |1 + _ H- - -j + 0(V<)| (52) 

for n = k + 1, and 

'-(V) - (-1)"" {V + + 0(V=)} (53) 

for n 7^ /c + 1. Hence, if the logarithmic dependence of (02n('Ti^)) on is 
neglected, the fc-th moment projects out the power corrections of dimension d = 
2k + 2: 

KU-i-^r^^^W^ + 0(V). (54) 

The most important effect of the scale dependence is to induce contributions from 
lower-dimensional operators. 

Using the above results we can derive explicit expressions for the power correc- 
tions to the moments, treating both V and l/m^ as small parameters. Neglecting 
terms of order {O2) /m?^, (04)/m^, V {OQ)/m^ and operators of dimension 

^This becomes explicit if these quantities are expressed in terms of scale-invariant conden- 
sates, see Appendix |c[ 
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larger or equal to eight, we find 

^2 



■' '0, power 



(J) 

power 



1 + v + 

{or 



1 



•' '2, power 



fe, power 



2V 

3 



V 



2 4 
k+1 



TC 

"6 



{0^"> 



{or'') 



V (or^)) 



(•/)\ 



+ ... 



+ . 



w - 1 w2 



2 



A: + 1 



(Of) 



+ 



Thus, the power corrections are parametrized by a set of fundamental nonper- 
turbative parameters V'^ {O^J)- 

We now apply this general formalism to calculate the leading power correc- 
tions in the SVZ approach The corresponding corrections {O^^ 
current correlators D^"^^ have been calculated, at next-to- leading order in the 
coupling constant, by several authors. They have been summarized and rewrit- 
ten in terms of scale-invariant condensates in Ref. [Q. For the corrections of 
dimension two and four, the results are P^-||17[] 



{or'') 
{or') 



-3 
27r2 



^ 13 «s 

3 TT 

11 a. 



a 



1 - 



18 TT 
TT 



TT 



UJ I 



+ 12 1 v; 



2 2 2 



{mulpulpu + ITT'jtpj^j) 

24 71 



— i'^k'ipM 



E 



4 



6 



a., 4 / 



^ (m^ + m^) + const. 



/24 TT 10 , 



{ml + m^) . 



(56) 



(Of) 
(Of) 

Here = as{ml) and = mqijn^). A summation over j = d,s is under- 
stood. The terms proportional to l/a^ arises from IR logarithms, which can be 
resummed into the running coupling constant by means of the RGE Note 
that (of) is defined only up to a constant (with respect to ml), which depends 
on the renormalization scheme. However, only derivatives of (Of ) appear in 
physical quantities. 

The most important power corrections of dimension six come from four-quark 
operators; the coefficient of the operator vanishes to leading order in ag 
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[^ , |i9| . Operators containing powers of the light quark masses are strongly 
suppressed and can safely be neglected. The coefficient functions for the four- 



quark operators were calculated to next-to-leading order in Ref. |50]. In the chiral 
limit there are no contributions to the J = correlator, and the contributions to 
the J = 1 correlator areQ 



647r^ ' 

+ 1677^ 

327r^ 



2L 



2An 



3 



1 + 
1 + 

1 - 



87 1 ■ 

V -L 

16 4 

1153 _ 95^ 
432 " 72 



TT 



{ol' + o 



ij ) 



TT 



TT 



^8 + O^*^) 



TT 



TT 



;i - 6L) 



TT 

2^^ 2 



E 
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+ Ojl + (Of + 



A8\ 



vr 



) 



(57) 



The dependence on resides in L = In / /i^ , and /i is an arbitrary renormal- 
ization scale. Since we work in the chiral limit, it is not necessary to specify the 
flavour labels; however, it is important that i ^ j- The four-quark operators, 
which are renormalized in the MS scheme, are defined as 



VI 



1] 



(5^ 



Since there exist no reliable estimates for these condensates, the traditional ap- 



proach is to apply the vacuum saturation (or factorization) approximation [|T5 



to relate them to the quark condensate. Unfortunately, this approximation is 
inconsistent with renormalization-group (RG) invariance |51|, 0. Moreover, phe- 



nomenological analyses indicate that the factorization approximation underesti- 
mates the contribution of four-quark operators. One introduces a fudge factor p 
to compensate for this. At leading order in as, one then obtains 



/fact 



2567r3 
27 



pasiip'ipy^ 



(59) 



The combination pas{ipip)'^ is treated as a scale- invariant phenomenological pa- 
rameter, whose value has been determined to be of order 2-6 x 10"^ GeV^ P-[pT[], 

^We have rewritten the result given in Ref. using Fierz identities. Note that, in the 
notation of Ref. pi, one has (Og"^^^) = Qtt^ 0(6). 
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corresponding to p ~ 3-8 at p = 1 GeV. The fact that phenomenology requires 
p ^ 1 is unsatisfactory and indicates large violations of the factorization hypoth- 
esis. 

An interesting and, to our knowledge, new alternative to the vacuum satu- 
ration approximation is to consider the large-/3o limit, which we have employed 
already in the calculation of the perturbative corrections. This limit respects 
exact RG invariance, yet reducing the dimension-six corrections to a single phe- 
nomeno logical parameter. To obtain it, we analyse expression {^7\j in the limit 
n/ — >■ oo. Note that there is dependence in the running coupling constant as 
well as in the flavour sums. We find 



(0+1), _ 647r4 fnfOsilJ^Y V8f,2 



{or'') = --f- [^''^) 0-(p^ 



TT 1 



+ — (1-6L) 



+ 0{l/nj) , 
(60) 

where the product njasifi"^) is formally of order n^. The two terms shown in 
parenthesis correspond to the diagrams depicted in Fig. ^. In the next step, we 
replace Uf by — f/^o and use the one- loop expression for the running coupling 
constant to find that 

(Or^^)ia.,e-,o = 128vr^ (l- ^ + ^) «'(/^') O^if^') , (61) 

where Ay = e^/^Aj^, and a(p^) has been defined in (^7]). The p dependence 
of a(p^) is exactly cancelled by the p dependence of the four-quark condensate 
. Note that the contribution (|6T1) vanishes in the factorization approximation. 
In this sense the above two approximations are orthogonal to each other. If we 
define 

where (ipip)'^ is evaluated at p = 1 GeV, we find that the results (0) and (|^) are 
of similar magnitude if e ~ 1-2%, corresponding to only a small deviation from 
vacuum saturation, which is certainly not excluded. Hence, if one accepts that 
the large-/?o limit gives a reasonable approximation to QCD, it provides a natural 
explanation of the empirical fact that the values of four-quark condensates are 
much larger than predicted in the vacuum saturation approximation. 

Starting from (|56| ) and using the /5-function and the anomalous dimension of 
the running quark masses from Appendix |C|, we find 

V(0r^)) = 6^|K,r(m^ + m,^), 



TC 



+ 6|K,f (m^ + mj). 
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Figure 4: Examples of leading contributions of four-quark con- 
densates IN THE LARGE-nj LIMIT. ThE FIRST DIAGRAM IS OF ORDER HjOs, 
THE SECOND OF ORDER (n/Os)^. 



V(Of) = 6|K,f K + mj). (63) 

The quantity (02°^^'') appearing in (|55|) is of order a1 and will thus be ne- 
glected. For completeness, we also quote the results for the derivatives of the 
dimension-six operator. They are 

(0+1), 21767r3 a, - ^ 

V (of +'^)large-/3o = 12871^ £ (V^V^) ^ . (64) 

To estimate the size of the various power corrections we use the standard 
values of the QCD parameters, which are collected in Appendix |C[ We find the 
following results for the nonperturbative parameters corresponding to J = + 1: 

^ ^ ' ^ -(1.42 ±0.27) X 10~^ 



ml 



mi 



(0.20 ±0.04) X 10-^ 
(9.2 ±5.6) X 10-^ 
(0.18 ±0.08) X 10-^ 

(3.3± 1.9) X 10'^ (65) 



The quoted errors reflect the uncertainty in the values of the quark masses and 
vacuum condensates. We note that the values of (04°^^'') and V (04°^"^'') are 
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dominated by the gluon condensate. For J = 0, we obtain: 




(0.88 ±0.17) X 10 



-3 



(Of) 



(0.10 ±0.02) X 10 



-3 




(1.92 ±0.72) X 10 



-3 



V(Of) 



0.01 X 10 



,-3 



(66) 




Most of the contributions are strongly suppressed by powers of the small quark 
masses and can safely be neglected. Contributions to the moments of order 10~^ 
are certainly not detectable given the uncertainties in the perturbative calcu- 
lation, and also given that any extraction of the moments will be affected by 
experimental uncertainties. Therefore, it is a safe approximation to consider the 
power corrections in the chiral limit, in which case only contributions with J = 1 
remain. The most striking prediction of the SVZ approach is that (Of ) = in 
the chiral limit, since there is no gauge- invariant operator of dimension two in 
QCD. As a consequence, the leading nonperturbative corrections are induced by 
the gluon condensate and are of order {0'^^)/m^ ~ 1%. We expect that correc- 
tions of dimension six are suppressed, relative to this, by about a factor 3, which 
is beyond the precision reachable in a realistic analysis. Therefore, we believe that 
the primary goals must be to measure ^^(m^) with a minimum contamination 
from power corrections, to test the absence of dimension-two operators in QCD, 
and to extract a value for the gluon condensate. We shall discuss some strategies 
how to pursue these goals in Sect. |^. Before, however, we will demonstrate how 
the inclusion of power corrections cures the IR renormalon problem. 

6 Cancellation of renormalon ambiguities 

The structure of power corrections simplifies greatly in the combined large-/5o and 
chiral limit. The nonperturbative parameters remaining in this limit are the ones 
needed to absorb the renormalon ambiguities in the perturbative calculations of 
Sects. I and |. We find 




V (Oj'))ia.ge-/3o = 1287r^£(V>V')'. 



(67) 
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All other contributions of dimension up to six vanish in this limit; in particular 
all contributions for J = 0. It is satisfying to observe that the potentially large 
power corrections survive in the large-/5o limit. From the numerical discussion of 
the preceding section, it follows that all other terms are of order 10~^ or less (in 
the appropriate units of rrir). 

Let us now demonstrate how the renormalon ambiguities of perturbation the- 
ory can be cured by adding nonperturbative corrections. To start with, consider 
the euclidean correlator In the combined large-/5o and chiral limit, we 

find from (^6|) and (|47|) 



D^'\Q') = l + ^ /-^^D(r)a(rQ^) + i^ + M^ + 0(Q-«). (68) 

The perturbative contribution is ambiguous because of the Landau pole in the 
running coupling constant a{TQ'^). The corresponding renormalon ambiguity has 
been given in (|3lD. To obtain a meaningful result, we have to require that, order 
by order in these ambiguities are cancelled by corresponding ambiguities 

in the definition of the nonperturbative parameters {O^^)- The philosophy is the 
following: The values of the nonperturbative parameters become meaningful only 
after one specifies a resummation prescription for the perturbative series, i.e. a 
prescription to regulate the Landau pole in the integral over r in (p^). Changing 
this prescription changes the values of the nonperturbative parameters in such a 
way that the total answer remains the same [Q. Hence, to order we have 
to require that 



Note that the momentum dependence of {0^\Q'^)), which is obtained from (pl|) 
by the replacement — > Q^, is precisely of the form required to fulfill this 
condition. It follows that 

a/^ gA = -1^ ^ -0.006 GeV^ , 

\ TT /ren VT^/^q 

- - 1 

el^^)^ = ~ 5 X 10~^ GeV'^ , (70) 

J ren 24:71^^^0 



A 



where we have used Ay = 0.461 GeV corresponding to as(m^) = 0.32. Numeri- 
cally, the ambiguity of the gluon condensate is about one third of its value; the 
ambiguity of the dimension-six condensate is about 3 x 10^^/e times its value, 
which amounts to 30% if e = 1%. 

Once the ambiguities cancel for the euclidean correlation function, they also 
cancel for the moments. We find that to order the conditions for this 
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cancellation read 

6 k — 2 
^ 3 

AAfr + ^<"°"y""° + -^'^'1"""° = o. (71) 

m° 3 m° 

where the perturbative ambiguities AMl*^^ have been given in (Q). These con- 
ditions are indeed satisfied with ([fO|). 



7 Discussion of the results 

We have presented a new analysis of the spectral moments measurable in hadronic 
r decays, with the main focus on the uncertainties inherent in the perturbative 
calculation of the leading terms in the OPE. Comparing fixed-order perturbation 
theory (which is known to order in this case) with a partial resummation of the 
perturbative series based on the analysis of so-called renormalon chains, we con- 
clude that the uncertainties in the calculation of the moments are of order a few 
per cent. We have also considered the resummation prescription of Le Diberder 
and Pich and find that it does not reduce this uncertainty in a significant way. 
These observations imply important limitations for tests of QCD at the level of 
power corrections. We have presented a model-independent description of such 
corrections and argued that they have to be included, for reasons of consistency, 
in order to cure certain ambiguities of perturbation theory, which are related to 
IR renormalons. For the power corrections of dimension four and six we have 
demonstrated in detail how these ambiguities can be absorbed into a redefini- 
tion of some nonperturbative parameters (vacuum condensates). Based on our 
results, we conclude that previous analyses of spectral moments in hadronic r 
decays |12|, 0, which were aiming at an extraction of power corrections up to 



dimension six and even eight, have underestimated the theoretical uncertainties. 
In the remainder of this section we shall reconsider the feasibility of extracting 
fundamental QCD parameters, such as the running coupling constant and some 
of the vacuum condensates, in the light of our results. 

Even in the ideal case in which the spectral moments Ai^^^ can be measured 
with arbitrary precision, and in which the perturbative contributions to these 
moments can be calculated to very high order, there is an unavoidable ambi- 
guity in the definition of power corrections, which results from the presence of 
renormalons in perturbative QCD. The renormalon ambiguities in the values of 
the gluon condensate and the dimension-six condensate have been given in (|TOD . 
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They amount to about 30% of the expected values of the condensates. How- 
ever, to achieve such a level of precision would not only require zero experimental 
errors, but also a control of the asymptotic behaviour of perturbation theory, 
which is lacking at present. One cannot trust the large-/3o limit to provide a very 
accurate description of this behaviour. To be conservative, we shall consider the 
difference between third-order truncated perturbation theory and the resummed 
results in the large-/9o limit as an estimate of the perturbative uncertainty. We 
will now discuss what this implies for measurements of as(m^) and some of the 
vacuum condensates. 

The most accurate way to determine the running coupling constant is to 
measure the total inclusive ratio defined in (||). We separate the perturbative 
contribution to this quantity from nonperturbative corrections by writing 

Rr = 3Sew (1 + ^port + S power ) • (72) 

The perturbative contribution, (5pert = 2Mo — 2M2 + M3 — 1, depends strongly 
on the value of the running coupling constant. In Fig. ^, we compare the exact 
order-a^ prediction for this quantity (Table |l]) to the result obtained from the 
resummation of renormalon chains (lower portion of Table ^). For comparison, 
we also show the result obtained using the resummation method of Le Diberder 
and Pich ^ (Table 0), which is routinely used in the determinations of from 
Rt-. The differences between these perturbative approximations are of order few 
to several per cent, increasing as the value of as increases. The quantity 

, _ ^ vx (or^)) ^ ^ V {or'') _ , {or'') 

Opower ~ \ ~^ o I ■? ~^ d fi 

V o / m° 

4v-Hv2)i^ + 2(l-V)^ + ...^-(1.4±0.5)F3) 
6 9 / mi 



contains the power corrections. In the SVZ approach, the largest contribution 
(about 50%) comes from dimension-six operators; the contribution of the gluon 
condensate is suppressed by two powers of as- The remaining corrections come 
mainly from quark-mass effects and from the quark condensate, each of which 
contribute about 25%. The fact that the power corrections to R^ are very small 
has been noted in Ref. 0. It is for this reason that measurements of R^ are 
believed to provide a reliable determination of the running coupling constant. 
Note that the uncertainty in the value of the power corrections is much less 
than the uncertainty in the perturbative calculation of 5port- Thus, it is the 
perturbative uncertainty which limits the accuracy in the extraction of as(m^). 
Even if one could determine 5pert with a very small error (in practice this error 
cannot be much smaller than 1% because of the theoretical uncertainty in (Jpower); 
there is a rather large uncertainty in the corresponding value of the running 
coupling constant, depending on which perturbative approximation one uses. 
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For S 



pert 

2^ r 



'.0.2, which is the value preferred by experiments |T6|, [T] 
0.303 from the resummation of renormalon chains, a. 



m. 



one obtains 
2^ ~ 0.337 



from the exact order-af calculation, and as{m'^) ~ 
of Le Diberder and Pich. Hence, we conclude that 



0.05 



0.353 from the resummation 



(74) 



is a reasonable estimate of the theoretical uncertainty. This is a factor 3 larger 
than the total error estimated in Ref. Q and quoted in the experimental 
analyses [11^, K7\ . 



0.4 



0.3 



CO 



0.2 



0.26 0.28 0.3 0.32 0.34 0.36 0.38 
as(m?) 

Figure 5: Different perturbative approximations for the quantity 

Spert- resummation IN THE LARGE-/?o LIMIT (SOLID LINE) , EXACT ORDER-af 
RESULT IN THE MS SCHEME (dASHED LINE) , RESUMMATION OF Le DiBERDER 
AND PiCH (dash-dotted LINE). 



Next consider the extraction of the leading nonperturbative parameter con- 
tributing to the moments, which is the gluon condensate or, more generally, the 
quantity (04'''''^^). From ([55| ) it follows that the moment Aif'^^^ is most sensitive 
to this parameter. Let us assume an ideal experiment, in which this moment 
can be extracted with high accuracy.^ Moreover, assume that as{m'^) has been 
determined independently from a measurement of Rr- As previously, we write 

Mf^'^^ = 1 + 77pert + r/power • (75) 

The theoretical predictions for the perturbative contribution, rjpert = Mi — 1, can 
be obtained directly from the numbers given in Tables |l|, |] and For values of 

^ In practice, this would require to subtract the contributions with J = 0. Moreover, 
what is accessible in experiments are not the spectral moments A^'^'' themselves, but rather 
the moments Rk of the r hadronic spectrum, which according to M) are combinations of the 
spectral moments. 
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asim"^) in the range from 0.30 to 0.34, the results differ by 1-3%. For the leading 
power corrections we obtain, neglecting terms of order 10""^, 

%ower = —1— + . . . ^ 0.2% - ^1-7^— (— C ). (76) 



3m4 V 18 TT / \ TT 

For {— G^) ~ 0.018 GeV^, the contribution of the gluon condensate is —1.1%, 
thus providing the dominant correction. The precision in the determination of 
the gluon condensate is limited by the theoretical uncertainty in the value of //pert- 
An uncertainty of 1% in the value of rjpert implies an uncertainty of 0.02 GeV^ 
in the value of the gluon condensate, which is of the same magnitude as the 
expected value of the condensate itself. We conclude that, unless one manages to 
suppress the perturbative uncertainty, the gluon condensate cannot be extracted 
with a precision of much better than 100%. 

Consider finally the extraction of power corrections of dimension six. An ideal 
observable to determine these would be the combination 

3A^(0+l) _ 2A<r'^ = 1 + Ppert + Ppower , (77) 

since it does not receive power corrections from {0^~^^^) and V (04''^^^). Neglect- 
ing a contribution of order 10~^ from V (02°^^''), we find 

Ppower = —2 7 h . . . . (78) 

In this case, the perturbative uncertainty in the value of ppert = 3M3 — 2M2 — 1 is 
of order 1-2%, leading to 6{0q^^^^) ~ 0.2-0.3 GeV^, which is larger than expected 
value of the condensate. Given the present level of control over the perturbative 
contributions to the moments, it is thus not possible to extract power corrections 
of dimension six or higher from hadronic r decays, even in an ideal experiment. 

To summarize, we believe that beyond extracting the running coupling con- 
stant asi'm'^) the goals of a realistic analysis of spectral moments in hadronic 
r decays should be (i) to test QCD, or more specifically the SVZ approach, by 
checking the absence of large dimension-two operators; (ii) to extract a value 
of the gluon condensate. To pursue these goals, we propose to measure differ- 
ences of the moments Rk defined in (^. The perturbative contributions to these 
differences start at order and are known exactly to order a^. Since these con- 
tributions are very small, the quantities we shall construct provide a more direct 
measurement of power corrections, which is very little affected by uncertainties 
in the value of ^^(m^). 

A quantity that is sensitive to the presence of non-standard power corrections 
of dimension two is 

2 

TTL 

1 . X 1 f / 10 s \ dRJs) 

g T 
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-Mnf -Mf + (79) 

9 o o 

Using (p!8D, we find that in the MS scheme the perturbative contribution to this 
difference is 

^pert = ^(^y + cJ^)\cJ^y + ...^ (2.8 ± 1.0)% , (80) 



80X77/ \ Tc J Vvr 



where 



241531 513 , 

^ 596310293 ^ im7 1935519 ^ 2565 

768000 1280 3200 ^ ^ 32 ^ ' v ; 

and we have used as = as{m^) = 0.32. Numerically, we find that every term in 
this series is of the same magnitude, and we have used the last term to estimate 
the uncertainty. If we perform a partial resummation of the series in the large-/3o 
limit and correct for the known pieces to order (this is a small corrections 
in this case), we find the larger value Dg"'^'^' = (5.5 ± 0.7)%, where the error is 
given by the ambiguity due to the nearest IR renormalon. The leading power 
corrections of dimension less or equal to six are 



2 J ml \3 18 / 

Z_^V^)<^+fl_l£v)<^ + .... (82) 

9 108 y m2 V 9 / m4 ^ ' 

In the SVZ approach the largest contribution comes from the gluon condensate, 
and we find Df^'"' ~ (0.8 ± 1.0)%. We conclude that 



(3.6 ± 1.3)% ; fixed-order perturbation theory, 
(6.3 ± 1.2)% ; resummed perturbation theory. 



^3) 



A non-standard dimension-two contribution of size (A/m^)^ ~ 3-8% (for A = 
300-500 MeV) and possibly negative sign could be large enough to change this 
value in a significant way. Such a contribution is forbidden in the standard SVZ 
approach, since there is no local, gauge-invariant operator of dimension two in 
QCD. We note, however, that the truncation of the perturbative series for D2 can 



fake such a term, because this series contains an UV renormalon at u = — 1 [43]. 
It would be interesting to see whether the resummation of renormalon chains, 
which resums this contribution in the large-/5o limit, improves the agreement 
with experiment. 
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Another interesting quantity, which can be used to measure the gluon con- 
densate, is 



10 35ew J TTT'i V4 mi J ds 

+ \Mf-^^M^H^MT-\Mf. (84) 



= _ _ d2 f-V + ■ ■ ■ - -(0-46 ± 0.10)% , (85) 



The perturbative contribution to this difference is 
with 

34527 243 , 

di = C 3 ~ 1.7436, 

6400 80 ^ 

77871081 2187 2 819369 , 1215 , 

d2 = TT^ C3+ (5^ 9.4508. 86 

^ 512000 2560 6400 64 ^ ^ ^ 

The resummation of renormalon chains gives a very similar result, Df°'^'^^ = 
— (0.55 ± 0.40)%. In both cases the perturbative contribution is very small, and 
the accuracy seems to be limited by the renormalon ambiguity. We also note that 
the coefficient of the UV renormalon at m = — 1 is very small in this case, so that 
one does not expect large uncertainties due to the truncation of the series. The 
leading power corrections to the quantity D4 are 



„ i3 v(or)) / 35 x(or^) 3(or^) 

^ AO ml \ 2A ) 2 m^^ 

/W_1Q7_ ,\ (Of) _ 7_ V(Of) 
V30 1800 / m2 24 

To good approximation, they are given by the gluon condensate. Combining 
them with the perturbative contribution, we obtain 

— 

Hence, a measurement of Di^ comes close to a "null measurement". If a value 
7^ can be established, it provides direct evidence for the existence of the 
gluon condensate. We would consider such a measurement more convincing than 
the existing fits |jT6|, |T^ of the gluon condensate from a moment analysis following 
the lines of Ref. |0 . 
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Appendices 

A Calculation of the integrals I^^^{mr) 

Le Diberder and Pich have introduced the integrals ||T2 



{nir) 



1 

27ri 
1 

2^ 



\x\=l 



X \ TC 



TT 



^9) 



The running couphng constant along the contour in the complex plane can be 
obtained from the solution of the RGE (|T^, which reads 



471 «,(0) 47r/?o ""aM levr^/^^ 



«,(¥^)-«,(0) +... , (90) 



with as{(p) = as{e"^ml). 

In the approximation where one uses the one-loop /3-function (setting jSi = 
/32 = . . . = 0), it is possible to perform the contour integrals explicitly. The result 
for n = 1 is 



where 



rriT 



4 

Wo 



TC 



g-(fc+i)/a ini Ei [(A; + 1) (1/a + in) 



" = ^"-^"^-^ = ln(mU^) ' 



H — arctan(7ra) > , (91) 



(92) 



and Ei(x) = y e* is the exponential integral. Note that the value of Im Ei(a;+ 
iy) is not affected by the pole at t = 0. The integrals with n > 1 can be obtained 
from the recursion relation 



/ A \ n-1 



1 



'n-D 



da 



-1 



n— 1 



k+l. 



m-r 



(93) 



Beyond this approximation, the integrals have to be performed numerically. 
However, in the case of the two-loop /3-function one can still obtain an approxi- 
mate analytic expression by writing 



1 I Pi as{^) 



asiif) a,(0) Airpo a,(0) 
77^ ■ 



167r2 



(94) 
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Table 10: Integrals /^+^ evaluated using the one-loop (upper portion), two- 
loop (middle portion), and three-loop (lower portion) /3-function. In the 
two-loop case, the approximate results obtained from (|9^) are given below 
the exact values. 





k = 


k = 1 


k = 2 


k = 3 


k = A 


k = 5 


(xlO- 




11.55 


9.73 


9.40 


9.24 


9.16 


9.10 


/|'+^ (xlO- 




12.08 


8.01 


7.62 


7.40 


7.28 


7.20 


1^+^ (xlO- 




11.67 


5.25 


5.22 


5.05 


4.95 


4.89 


1^+^ (xlO- 




11.44 


9.17 


8.91 


8.75 


8.67 


8.61 


approx. 




11.57 


9.37 


9.04 


8.88 


8.78 


8.72 


1^+^ (xlO- 




11.54 


6.68 


6.65 


6.44 


6.34 


6.28 


approx. 




11.79 


6.95 


6.74 


6.54 


6.43 


6.36 


Jf+^ (xlO- 




10.56 


3.34 


4.12 


3.92 


3.86 


3.82 


approx. 




10.86 


3.51 


4.03 


3.89 


3.82 


3.78 


(xlO- 




11.42 


9.11 


8.87 


8.71 


8.63 


8.57 


1^+^ (xlO- 




11.48 


6.58 


6.59 


6.38 


6.28 


6.22 


/|+^ (xlO- 




10.46 


3.21 


4.08 


3.85 


3.80 


3.76 



This leads to the simple relation 

where in the one-loop integrals one uses 

l + ^a)a (96) 

instead of a. This provides a very good approximation to the exact numerical 
results. 

In Table [1^, we compare the various approximations for some of these integrals 
in the case of ^^(m^) = 0.32 (in the MS scheme). 

B Asymptotic behaviour of the functions Wk(T) 

In order to perform numerical integrations with the distribution functions Wkir) 
defined in (|35D -(^7|), it is convenient to use their asymptotic behaviour for large 
and small values of r. For r ^ 1, we find 

, . /, , X , , 5fc + 16 Inr \ 1 

Wk(r) = — A; + 1 . + 



36(fc + 2)2 6(A; + 2)/r 
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^ 7fc + 33 , Inr \ 1 ^ ^/_-3n 



which is vahd for all A; > 0. For r <^ 1, on the other hand, we obtain 

Woir) = f { (4 - 3C(3)) r - ^ + (1 - i Inr) + 0(r^)} , 
W,ir) = f I (5 - 6C(3) - ^ Inr) + - Inr) r^ + 0(r^)| , 

l^,(r) = f - - i + i 1-'-) + 0(r^)} , (98) 

as well as 

for k > 3. 



The asymptotic behaviour of the distribution function wj:,{t) in (|29|) can be 



derived from the above results by taking the limit k ^ 00. For r ^ 1, we find 

(.nr + 5)i-5(.„r + ^)l + 0(.-), (100) 



Wo(r) = 8r2 + 8 ( - In r - 1 ) r^ + O(r^) . (101) 



whereas for r -C 1 

wr^ir) = 8r2 + 8 

.3 

C QCD parameters 

For the running quark masses at the scale /iq = 1 GeV (in the MS scheme), we 
use 

m„(/io) = 5± 1 MeV, 
mdifJ^l) = 9 ± 1 MeV , 

rusif^l) = 178 ± 18 MeV . (102) 



The value for has been obtained in a recent QCD sum rule analysis . The 
scale dependence of the running masses is governed by the RGE 

^^^=7m[asi^^')]mli^^'). (103) 

At two-loop order, 

7n^(«.) =70^ + 71 (^)' + .. - , (104) 
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with coefficients 70 = —8 and 71 = —364/3 (for Uf = 3) [^. At tliis order, tlie 
solution of the RGE is 



m.il^ ) = m^i^^,) f ^ j (1 + ^^^ ^ + • " " 



(105) 



This can be used to calculate mq(m^). In the numerical analysis we use ^^(m^) = 
0.32, unless otherwise specified. 

In the calculation of the power corrections in Sect. |^ we also need the elements 
Vuj of the Cabibbo-Kobayashi-Maskawa matrix. We neglect Vut, so that | V^dp + 
I Ktsp = 1, and use 



I Kd|= 0.9753, I Ksl = 0.2209. (106) 

Finally, we give the definition of the scale-invariant condensates of dimension 
four, which are relevant to our analysis. In the MS scheme, and to next-to-leading 



order in the coupling constant, they are defined as ^6] , ^ 
^ = f 1 + ^) ^ (G,.G-(/.^)) 



7r / \ 9 TT / TT 

16 asifi"^) A , 19 a^(/i2 



,1 + ^^^ E^fc(/^')(V^^^^(/^')) 

1 ^1 



37r^ V 3 77 . ^ 
{m,^P,^P,) = m^ifi') {^Mf^')) + ^ - §) • (107) 

For the scale-invariant quark condensates we write {nii'ipi'ipi) = —mil^l) Hi 
take the standard values 11151 



= Kd = (230 ± 30 MeV)^ , = 0.65k„ . (108) 

For the scale- invariant gluon condensate and for the four-quark condensate we use 
an average over some recent determinations |0-||lT[, however with conservative 
errors: 

= (0.018 ±0.009) GeV\ 

TT / 

pasitpt/j)'^ = (3.5 ± 2.0) GeV*^ . (109) 
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